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Целью работы является анализ притягивающих множеств, фор- 
мируемых при фрезеровании концевыми фрезами. Предвари- 
тельно рассмотрено формирование стационарных траекторий 
периодических деформационных смещений инструмента отно- 
сительно заготовки и их устойчивости. В данном случае движе- 
ния образуют притягивающее множество, задаваемое траекто- 
риями периодических изменений толщины и ширины срезаемо- 
го слоя каждым зубом фрезы с учетом деформационных смеще- 
ний. В отличие от ранее рассмотренных случаев, в статье анали- 
зируются притягивающие множества, формируемые при потере 
устойчивости стационарной траектории. Показано, что в зави- 
симости от параметров системы и технологических режимов в 
динамической системе фрезерования могут формироваться при- 
тягивающие множества типа предельных циклов, инвариантных 
торов и странных (хаотических) аттракторов. При этом анали- 
зируются два случая. В первом параметры системы (прежде 
всего — скорости резания) позволяют пренебречь изменениями 
коэффициентов в дифференциальных уравнениях в пределах 
импульсной реакции системы. Во втором случае параметры 
системы изменяются в пределах импульсной реакции системы 
— ив ней образуется дополнительный источник параметриче- 
ского самовозбуждения. Значительное внимание в статье уделе- 
но анализу бифуркаций притягивающих множеств в параметри- 
ческом пространстве: дается общий анализ и приводятся приме- 
ры. Притягивающие множества анализируются под углом зре- 
ния их влияния на параметры качества изготовления деталей. 


Ключевые слова: процесс фрезерования концевыми фрезами, 
притягивающие множества, бифуркации, управление, парамет- 
ры качества изготовления деталей. 
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Введение. В последние два десятилетия при изучении свойств эволюции и самоорганизации широко исполь- 


зуется синергетическая парадигма [1-4]. Ее применение для управления сложными нелинейными объектами нашло 


свое отражение в работах [5-8]. При создании систем управления динамической системой резания также используют- 


ся основы синергетической теории управления [7-14]. В этом случае управление, в том числе на основе построения 


программы ЧПУ, включает определение желаемых траекторий формообразующих движений и соответствующих им 
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траекторий движения исполнительных элементов станка. При решении этой задачи принципиально важно знать свой- 
ства тех притягивающих множеств, которые самостоятельно образуются в окрестности формообразующих движений 
инструмента относительно обрабатываемой детали. Необходимо также уметь управлять этими множествами, влияю- 
щими на качество поверхности, формируемой при резании. Наконец, раскрытие свойств притягивающих множеств, 
формирующих сигнал виброакустической эмиссии, открывает новые пути построения информационных моделей для 
диагностирования процесса [15]. 

Изучение притягивающих множеств связано с рассмотрением динамической системы резания, изучению ко- 
торой посвящено множество известных исследований [16-32]. Фрезерование является наиболее сложным процессом 
резания [33—46]. Это обусловлено его нестационарностью. Параметры длины и толщины срезаемого слоя каждым ре- 
жущим лезвием периодически изменяются. Поэтому в подвижной системе координат, движение которой определяется 
траекториями исполнительных элементов, стационарным, установившимся состоянием является не точка равновесия 
(как при точении), а замкнутая траектория. В связи с этим уравнение динамики имеет периодически изменяющиеся 
коэффициенты. 

Рассмотрим процесс фрезерования на станках, имеющих до пяти координат управления при обработке дета- 
лей, параметры жесткости которых существенно изменяются вдоль траектории движения инструмента [44—46]. В 
настоящей статье, в отличие от известных исследований, учитываются несколько источников самовозбуждения и свя- 
зи, обусловленные взаимодействиями передней и задней поверхностей инструмента с деталью. Здесь можно рассмат- 
ривать два случая. Первый — обработка с малыми частотами вращения шпинделя, когда в системе в пределах им- 
пульсной реакции параметры можно считать замороженными. Второй — обработка с большими частотами вращения, 
когда параметры нельзя считать замороженными. 

Математическая модель. Рассмотрим попутное фрезерование боковыми лезвиями прямозубой концевой 
фрезой недеформируемой заготовки (рис.1). 
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Рис. 1. Схема взаимодействия подсистем инструмента и обрабатываемой заготовки: деформации при движении стола 
в направлении образующей Г. (а), (6); деформационные смещения в подвижной системе координат (6), (г); силы резания 
в координатах состояния (д); учет матриц инерционных, диссипативных и упругих коэффициентов подсистемы инструмента (е) 


Деформации рассматриваются в плоскости, нормальной к оси инструмента. Считается, что ширина В обраба- 
тываемой заготовки постоянна, то есть В = соп5!. Толщина срезаемого слоя определяется а(ф) = а, (ф)-Х\ ($). Здесь 
а, (ф) — кинематическая составляющая толщины, зависящая от угла ф поворота фрезы. Она не учитывает упругие 


деформационные смещения инструмента и является периодической функцией угла. 
Введем в рассмотрение понятие — стационарной — траектории  деформационных смещений 


Х" ($) ={Х! (ф),Х. ($) }. Она характеризует установившиеся упругие деформационные смещения инструмента и яв- 


ляется периодической функцией ф с тем же периодом. Если в системе сформировалась стационарная траектория, то 


для нее справедливо: а(ф) =а(ф-1{Аф), а, (ф)=а,($ф-1ф), Х”(ф=Х"(ф-146),1= 0,1,2,... Рассматривается случай 
ГА 
О = сопя, где О — циклическая частота вращения фрезы. Тогда ф= [04 = От. Следовательно, справедливо: 
0 
а( =а((—Т), а. =а,.@&=—Т), х@Ф=хХ“аб=т),Е=0,1,2... 
Рассмотрим две системы координат. Первая имеет начало координат в центре фрезы без учета упругих де- 


формационных смещений. В этой системе Х (0 ое _ =) рассматриваются деформации центра фрезы. В 


частности, при продольном фрезеровании система движется по направлению Г, со скоростью у (рис. 1). В общем 
случае траектория центра вращения фрезы без учета ее упругих деформационных смещений определяется двумя фа- 
зовыми траекториями (УФ = ав /@,[) и п _ = АГ, / 4, Г.). Это программа ЧПУ станка. 

Вторая система привязана к кинематической траектории движения вершины режущего лезвия. В данной си- 


Т 
стеме рассматриваются деформационные смещения вершины режущего лезвия Х ={Х,,Х,} с \®. Связь между М и 


АЗ определяется 


хи — 5х ра 5х9 (1) 
где 
$1И с0$ ЗП —с0$ 
р (2) 
—с0$ф зШф с0зф  зШФф 


Так как <> = соп5т, то операторы 5(61) и 5 (ОГ) есть периодические функции времени. При рассмотрении 
системы в пространстве № фрезерование имеет свойства, аналогичные точению, но обработка ведется с периодически 


изменяющейся толщиной срезаемого слоя. Пусть в этом пространстве заданы силы, действующие на режущее лезвие. 
Тогда в пространстве № 


2 
пни чех = РР (3) 
ГА ГА 





т 
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т 0 п И. | Чл С21 
где т= ‚й= = — матрицы инерционных, диссипативных и упругих коэффициен- 


С = 
) 
О т По № С12  С22 
о > о Т 

тов инструмента; Е” (1) ={Е (6, Е (6) 

Параметры т, й и с можно определить по правилам, изложенным, например, в [7, 8, 47—49]. 

При полной симметрии инструмента матрицы й, с не только симметричные, но и диагональные [50]. Кроме 

1 

того, они имеют равные коэффициенты и не изменяются при переходе от пространства М к №. Использование си- 
нергетического подхода требует выяснения законов взаимодействия системы со средой [3-7] (в нашем случае — с 
процессом резания). Поэтому необходимо силы резания представить в координатах состояния. Они определяются 
суммой сил Е” =Е® + Е (рис. 1). 

Следуя [7, 8, 47—49], можно утверждать: 
— силы возрастают при увеличении площади срезаемого слоя; 
— имеет место запаздывание между изменениями площади и вариациями сил; 


— при увеличении скорости резания наблюдается уменьшение сил, которое в скоростном диапазоне (40-180) м/мин 


можно считать экспоненЦиальным. 


я, 
Сила Е“ резко увеличивается при увеличении скорости движения инструмента в сторону заготовки. Зави- 


симость этой силы от колебательной скорости характеризует нелинейное демпфирование. Кроме этого используется 


гипотеза о неизменности ориентации сил в установившемся состоянии. Она определяется для силы Е и УГЛОВЫМИ 


1 1 ь 
коэффициентами %х! и Х›. Для силы ЕЁ ии х' и хо . Этим свойствам соответствуют следующие представления сил 


ЕЕ, ЕО}. 


ах 
Е (Хи, Хо) = Хиро + пехр(-о (У — | Ва, (®-У1+хФЕ® -ехр[-а5 АХ, / а! 


ах 
Е) (Хи, Хо) = Хоро + пехр(-о (У — | В -[а,(®)-У,]+хЕ® -ехр[ сах, / ай; 


(4) 
АУ 
П-А+И=Х,; 
4 
АУ 
ВХ, 
4 
где Г, Т5 —Щ постоянные времени, моделирующие запаздывание сил по отношению к смещениям, с, причем Т,›Т 


(таким образом учитывается фазовый сдвиг между вариациями сил в направлениях Х\,Х,); ро — давление стружки 


2 0 & 
на переднюю грань инструмента, кг/мм’; Е ) — сила, действующая на заднюю грань в состоянии равновесия, кг; 
о, &> — параметры, определяющие крутизну изменения сил, с/мм; У — скорость резания. 


Модель (3), дополненная (4), описывает динамику фрезерования. 

Здесь обработка ведется без перекрытия контактов зубьями фрезы. 

Бифуркации динамической системы при отсутствии параметрического самовозбуждения. При малой 
скорости резания изменениями параметров системы (3) с учетом (4) в пределах импульсной реакции можно прене- 


бречь и рассматривать систему с замороженными коэффициентами [51]. В этом случае необходимо анализировать 
Й 0 а“) [52—59 Р 
свойства системы для а, = (а; ‚а, ’) [ ], то есть при постоянных параметрах на этом отрезке. Рассматриваемая 


система моделирует два взаимосвязанных колебательных контура, которые имеют два источника самовозбуждения. 
Один обусловлен запаздыванием сил, второй вызван существованием участка, в котором увеличению скорости соот- 
ветствует уменьшение сил. Поэтому в зависимости от параметров могут существовать три притягивающих множе- 
ства: точка равновесия, предельный цикл, двумерный инвариантный тор. 


Рассмотрим систему, параметры которой приведены в табл. 1. 
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Таблица 1 
Параметры динамической модели инструмента 


2 
т, кг.с/мм й, кг.с/мм с, кг/ММ 





500 0,0 
0,0 500 


Рассмотрим процесс продольного попутного фрезерования стали 20Х. Основные геометрические параметры 


0 0 ы ы х 
инструмента: у =6 , а=20’. Частота собственных колебаний инструмента без динамической связи, формируемой 
процессом резания: © = 710,0 Гу. Частоты системы с учетом динамической связи, формируемой процессом резания: 

(2) _ (2) _ 
«р =1540 Гу, > =840 Гц. 
Образование двух частот обусловлено тем, что за счет реакции со стороны резания, матрицы упругости ста- 
новятся, во-первых, несимметричными, во- вторых, не диагональными. 
2 
Давление стружки на переднюю поверхность инструмента — ру = (45,0 —60,0) кг / мм”. 


Коэффициенты ориентации силы: )\ =0,6, Х› =0,8, хо = хо =0. Они удовлетворяют условиям норми- 


7. 2 ь 
ровки — например, (%/)“ +(%5)` =1. Чтобы варьировать основные параметры динамической характеристики процесса 


резания, изменялись технологические режимы: 
— скорость резания в пределах 30—110 м/мин; 
— ширина срезаемого слоя в пределах 5,0—20,0 мм. 

Величина подачи на зуб неизменна и равна 0,1 мм/об. Проанализируем диаграмму (рис. 2) в плоскости пара- 
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Рис. 2. Бифуркационная диаграмма в плоскости (1 ,р): А — О — О — область устойчивости равновесия; 





В — О — С — область формирования двумерного инвариантного тора; А — О — В — область формирования 
орбитально-асимптотически устойчивого предельного цикла на низкой частоте; С — О — ОР — область формирования 
орбитально-асимптотически устойчивого предельного цикла на высокой частоте 








Здесь принято Т = 0,5Т,. Для точек 1, 2, 3, 4 на рис. 3 приведены фазовые траектории в сечениях фазового 


пространства. 
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Рис. 3. Фазовые траектории системы, соответствующие точкам 1 (а), (6); 2 (в), (г); 3 (9), (е); 4 (ж), (з) на рис. 2 


Все фазовые траектории сверху вниз соответствуют точкам 1, 2, 3, 4 на рис. 2. Кроме того, на рис. 2 приведе- 
ны также области О-разбиения для линеаризованного в окрестности равновесия уравнения. Они показаны пунктир- 
ными кривыми. На участках А — О — О пунктирные кривые совпадают с границами бифуркационной диаграммы, 
поэтому они сливаются с границами областей О- разбиения. Они определяют четыре области: 

—щ область, в которой все корни характеристического полинома линеаризованной системы расположены в отрица- 
тельной комплексной полуплоскости А — О — О; 
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— области, в которых пара комплексно сопряженных корней расположена в правой комплексной полуплоскости А — 
О Ежи Е— О— В; 
— область, в которой две пары комплексно сопряженных корней расположены в правой комплексной полуплоскости 
Е О— Е( рис. 2). 

Как видно (рис. 2), границы областей О-разбиения не совпадают с границами бифуркаций притягивающих 
множеств в параметрическом пространстве. Для объяснения этого можно рассмотреть фазовые траектории (рис. 3), 
соответствующие точкам 1, 2, 3, 4. Временные диаграммы для точки 2 приведены на рис. 4. 
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Рис. 4. Пример преобразования двухчастотного процесса в одночастотный, 
то есть в орбитально-асимптотически устойчивый предельный цикл 


Здесь видно, что двухчастотный процесс со временем преобразуется в одночастотный, то есть имеет место затя- 
гивание высокочастотных колебаний к низкочастотным. Аналогичная ситуация возникает и в области С — О — Е. 
Однако здесь низкочастотный процесс стягивается к высокочастотному. Эти эффекты связаны с нелинейной диссипа- 
цией, которая учитывается дополнительной связью, формируемой в области контакта задней грани инструмента с за- 
готовкой. Они определяются силой ЕР (2) (М /ай=Е (0) ‚ехр[-оаАХ\ / а (хо, хо }. Приведем также спектральные 
характеристики колебаний (рис. 5) в направлении Х/ для точек 2, 3, 4 (левые иллюстрации), а также траектории уста- 
новившихся колебаний в плоскости Х\ —Х. (правые иллюстрации). 
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Рис. 5. Изменение спектров установившихся колебаний и их ориентация в плоскости Х| -Х> 


для точек 2 (а), (6); 3 (в), (г); 4 (д), (е) на рис. 2 


Анализ влияния упруго-диссипативных параметров подсистемы инструмента на формируемые многообразия 
показывает: качественно бифуркационная диаграмма соответствует рис. 2. Варьирование с и параметров динамиче- 
ской характеристики процесса приводит к смещениям корней характеристического полинома в комплексной плоско- 
сти. При этом корни могут сближаться. В этом случае область, в которой формируется двумерный инвариантный тор, 
вырождается. Таким образом, при изменении параметров упругости и диссипации подсистемы инструмента происхо- 
дят изменения притягивающих множеств системы. В наших исследованиях образование хаотических аттракторов для 
рассматриваемой системы не обнаружено. Типичными притягивающими множествами в системе являются предель- 
ные циклы, в редких случаях — инвариантные торы. 

Ситуация принципиально меняется, если в системе дополнительно существует параметрическое самовозбуж- 
дение. 


Влияние параметрического самовозбуждения. По мере увеличения частоты вращения инструмента насту- 
пает момент, когда время периодического изменения коэффициентов становится соизмеримым с импульсной реакци- 
ей системы. На практике такие условия соответствуют скоростному фрезерованию. К аналогичному эффекту приво- 
дит и уменьшение изгибной жесткости инструмента. В этом случае в системе (3), (4) дополнительно образуется ис- 
точник самовозбуждения, обусловленного периодическими изменениями параметров в динамической связи. Рассмот- 
рим систему, параметры которой без периодической составляющей соответствуют точке 1 на рис. 2. В этом случае без 
периодического изменения параметров в окрестности единственного равновесия системы формируется притягиваю- 
щее множество типа предельного цикла на частоте, близкой к © =300 Гц. За счет параметрического самовозбужде- 
ния образуется дополнительный механизм потери устойчивости. Тогда в колебаниях становится заметным взаимодей- 


ствие параметрических колебаний и автоколебаний. 
Моделирование стационарных траекторий, иллюстрирующее это взаимодействие, приведено на рис. 6. 
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Рис. 6. Притягивающие множества в окрестности равновесия по мере увеличения частоты периодического изменения 


параметров: <> =150 Ги (а); @ =180 Ги (6); @= 210 Гц (в); Ф=230 [ци (г); О=260 Ги (9); 
(е); =300 [ц (ж); Ф=320 Ги (3) 


На данной иллюстрации показано изменение колебательных смещений инструмента в плоскости, нормальной 
к оси инструмента. Как видно, по мере увеличения частоты изменения параметров наблюдается синхронизация, в том 
числе дробного порядка, автоколебаний с периодическими движениями, генерируемыми переменностью параметров в 
уравнении (3) с учетом (4). 

При определенных частотах, зависящих от отношения частоты периодического изменения параметров и соб- 
ственных частот линеаризованной системы, образуются стационарные траектории типа инвариантного тора ( рис. 6, 
а, ж). 

Образуется также стационарная траектория, которая фактически состоит из двух чередующихся торов 
(рис. 6, г). 

При переходе от одного детерминированного стационарного состояния к другому наблюдается вначале раз- 
мытие траекторий ( рис. 6, 6, д, з), а затем образование хаотических колебаний (рис. 6, в, е). 

Характерно, что переход к хаотическим колебаниям осуществляется через бифуркации удвоения периода по 
сценарию М. Фейгенбаума [60]. Для этого проанализируем временные и спектральные характеристики движений ин- 


струмента (рис. 7). 
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Рис. 7. Временные и спектральные характеристики колебательной скорости в направлении Ху, 


соответствующие рис. 6 (а), (6), (в) 


На приведенных иллюстрациях видны изменения спектральных свойств колебаний и образование через удво- 
ение периода хаотических аттракторов. 

Обсуждение результатов. Фрезерование концевыми фрезами, как правило, сопровождается образованием 
различных притягивающих множеств, влияющих на процесс обработки и формируемую резанием поверхность. Как 
показано выше, в зависимости от параметров взаимодействующих подсистем, технологических режимов и геометрии 
инструмента возможно образование притягивающих множеств типа предельных циклов, инвариантных торов и хао- 
тических аттракторов. Поэтому при управлении процессом фрезерования необходимо не только обеспечивать требуе- 
мые траектории движения исполнительных элементов станка, но и заданные свойства системы, которые характеризу- 
ются типом притягивающего множества. 

Известно, что, например, интенсивность изнашивания инструмента зависит от параметров его автоколебаний 
[7, 8]. Имеются данные о том, что возбуждение колебаний в зоне резания в отдельных случаях вызывает уменьшение 
интенсивности изнашивания инструмента, в других — приводит к повышенному износу. Известно также, что интен- 
сивность изнашивания зависит от мощности необратимых преобразований подводимой от станка энергии в зоне обра- 
ботки. Повышение мощности необратимых преобразований всегда увеличивает интенсивность изнашивания. 

Если рассматривать динамическую структуру процесса фрезерования, то необходимо выделить те силы, кото- 
рые характеризуют необратимые преобразования. Во-первых, это диссипативные силы, которые формируются сим- 
метричной составляющей матрицы скоростных коэффициентов. Кососимметричные составляющие, характеризующие 
гироскопические силы, на виртуальных перемещениях работу не совершают. Поэтому в матрице скоростных коэффи- 
циентов работу совершают силы, формируемые симметричной составляющей матрицы скоростных коэффициентов. 
Кроме этого в матрице упругости работу совершают циркуляционные силы, которые определяются кососимметрич- 
ной составляющей матрицы упругости. Таким образом, при анализе влияния колебаний на интенсивность изнашива- 
ния необходимо рассматривать сумму сил, формируемых симметричной составляющей матрицы скоростных коэффи- 
циентов и кососимметричной составляющей матрицы упругости. Именно перераспределение этих сил в нелинейной 
системе может характеризовать уменьшение мощности необратимых преобразований в зоне резания. Этим, на наш 
взгляд, объясняется существование периодических траекторий движения инструмента, для которых возможно умень- 
шение суммарного значение мощности при определенных амплитудах и пространственной ориентации колебаний. 
Следует отметить, что в зависимости от частоты имеет место оптимальное значение автоколебаний, при которых 
стойкость максимальна. 

Колебания оказывают противоречивое влияние на формируемую при резании поверхность. Здесь принципи- 
альное значение имеет ориентация колебаний в пространстве. На рис. 8 приведен пример изменения ориентации коле- 


баний для рассмотренной выше системы в зависимости от частоты. 
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Рис. 8. Изменение ориентации колебаний инструмента в плоскости в зависимости от частоты возбуждения колебаний 


Здесь в зависимости от ориентации автоколебаний в пространстве они могут приводить к уменьшению вели- 
чины шероховатости формируемой поверхности. Колебания ограниченной амплитуды в направлении образующей 
формируемой поверхности могут существенно улучшать показатели качества формируемой поверхности. Образова- 
ние автоколебаний в направлении образующей всегда ухудшают шероховатость поверхности. 

Наконец, все нелинейные функции в математической модели не являются симметричными в вариациях отно- 
сительно точки равновесия. Поэтому колебания вызывают смещение точки равновесия системы за счет образования 
динамического смещения. Это один из механизмов влияния колебаний на геометрический размер формируемой дета- 
ли. 

Заключение. При управлении процессами обработки на станках необходимо решать две проблемы. Первая 
связана с проектированием траекторий движения исполнительных элементов станков и обеспечением их асимптоти- 
ческой устойчивости. Вторая — с обеспечением заданных свойств траекторий, которые характеризуются типом и па- 
раметрами притягивающих множеств, формируемых в окрестностях этих траекторий. Для управления типом форми- 
руемых притягивающих множеств необходимо не только изменять технологические режимы, но и выбирать соответ- 
ствующие параметры динамической структуры процесса — например, геометрические параметры инструментов. 
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